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The Quasi-Elastic (QE) contribution of the nuclear inclusive electron scattering model developed 
in Q is extended to the study of electroweak Charged Current (CC) induced nuclear reactions, at 
intermediate energies of interest for future neutrino oscillation experiments. The model accounts 
for long range nuclear (RPA) correlations, Final State Interaction (ESI) and Coulomb corrections. 
RPA correlations are shown to play a crucial role in the whole range of neutrino energies, up to 500 
MeV, studied in this work. Predictions for inclusive muon capture for different nuclei, and for the 
reactions pr)X and ^^C{vci,e~)X near threshold are also given. 


I. INTRODUCTION 

Neutrino properties have been object of much interest as long as they could provide hints of physics beyond the 
standard model. A sensitive way to study the mass of the neutrino is by means of neutrino oscillations. One of 
the various experiments devoted to this topic is the atmospheric neutrinos detection carried out in Kamiokande and 
Superkamionde, which data have given evidence on > Vr oscillation with 10“^ < ;< 10“^ and almost maximal 

mixing angle, see [3| for a review. Once this phenomena have been firmly stablished new questions arise, such as the 
role of three flavour oscillations and the precise determination of the values of neutrino masses and mixing parameters 
Q- For the obtention of accurate results in these new experiments it is necessary to keep under control the sources 
of systematic error. Two of the major sources of systematic errors in the sub-GeV samples of SK experiments are the 
charged and neutral-current cross sections 0 . Thus, if we want to cope with the requirements of this new experiments 
precise nuclear interaction models must be used. 

Any model aiming at describing the interaction of neutrinos with nuclei should be firstly tested against the existing 
data on the interaction of real and virtual photons with nuclei. At intermediate energies (nuclear excitation energies 
ranging from about 100 MeV to 500 or 600 MeV) three different contributions should be taken into account: i) 
Quasi-Elastic (QE) processes, ii) pion production and two body processes from the QE region to that beyond the 
A(1232) resonance peak, and iii) double pion production and higher nucleon resonance degrees of freedom induced 
processes. The model developed in P| (inclusive electro-nuclear reactions) and 0 (inclusive photo-nuclear reactions) 
has been successfully compared with data at intermediate energies and it systematically includes the three type of 
contributions mentioned above. The building blocks of this model are: i) a gauge invariant model for the interaction of 
real and virtual photons with nucleons, mesons and nucleon resonances with parameters determined from the vacuum 
data, and ii) a microscopic treatment of nuclear effects, including long and short range nuclear correlations ESI, 
explicit meson and A(1232) degrees of freedom, two and three nucleon absorption channels, etc. Finite size effects 
are computed from a Local Fermi Gas (LEG) picture of the nucleus, which is an accurate approximation to deal with 
inclusive processes which explore the whole nuclear volume The parameters of the model are completely fixed 
from previous hadron-nucleus studies: pionic atoms, elastic and inelastic pion-nucleus reaction^ A—hypernuclei, etc. 
[3- The photon coupling constants are also determined in the vacuum. Thus the model of Q and Q has no free 
parameters, and hence these results are predictions deduced from the nuclear microscopic framework developed in 
and 0. In this talk, we show an extension of the nuclear inclusive QE electron scattering model of 0, including the 
axial CC current, to describe neutrino and antineutrino induced nuclear reactions in the QE region. We will not show 
here many details of the model, for a detailed discussion we refer the reader to 0. 


II. INCLUSIVE CROSS SECTION 

We will present here the general formalism focusing on the neutrino Charged-Current (CC) reaction 

ui{k) + Az ^ I-{k') + X (1) 

though the generalization of the obtained expressions to antineutrino induced reactions or inclusive muon capture in 
nuclei is straightforward. 
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FIG. 1: Diagrammatic representation of the neutrino self-energy in nuclear matter. 


The double differential cross section, with respect to the outgoing lepton kineniatical variables, for the process of 
Eq. is given in the Laboratory (LAB) frame by 


dn{k’)dE[ |fc I 47r2 


( 2 ) 


with L and W the leptonic and hadronic tensors, respectively. 

The hadronic tensor includes non-leptonic vertices and corresponds to the charged electroweak transitions 

of the target nucleus to all possible final states; it is completely determined by six independent, real Lorentz scalar 
structure functions Wi{q^) , i = 1,..., 6. 

We follow here the formalism of P|, and evaluate the self-energy p) of a neutrino with helicity r inside of a 

nuclear medium of density p. Diagrammatically this is depicted in Fig. ^ After summing over helicities, we get 


E^{k]p) 


8iG 

V2M^ 


d^q {q; p) 

(27r)4 A:'2 — mf + ie 


( 3 ) 


where n(^(g) is the W"*'—boson self-energy in the nuclear medium. 

The neutrino disappears from the elastic flux, by inducing one particle - one hole (Iplh), 2p2h ... excitations, 
A(1232)—hole (Ah) excitations, or creating pions, etc... at a rate given by 


T{k;p) 


-pImE^(fc;p) 


( 4 ) 


To evaluate the imaginary part of Ei, we use the Cutkosky’s rules, and we cut with a straight vertical line (see Fig. 
^ the intermediate lepton state and those produced by the IT—boson polarization (shaded region). Those states are 
then placed on shell by taking the imaginary part of the propagator, self-energy, etc. Thus, we obtain 


Iml]i,(fc) 


8G 

V2M^ 


d^k' 0(gO) 
(2^)3 2E[ 




( 5 ) 


for > 0. 

Since TdtdS provides a probability times a differential of area, which is a contribution to the (yi^l~) cross section, 
we have 


da = r{k\p)dtdS 


— ■^ImE^(fc; p)dtdS 

rC 


— j^ImEi^(fc; p)d^r 


( 6 ) 


so the nuclear cross section is given by 


1 

W\ 


ImEi^(fc; p{r))d^r 


a = — 


( 7 ) 
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FIG. 2: Diagrams of some processes contributing to the self-energy. 


where we are considering a function of the nuclear density p(r) at each point of the nucleus and we integrate 
over the whole nuclear volume. We assume LDA, which, as shown in is an excellent approximation for volume 
processes like the one studied here. Coming back to Eq. o we can compare it with Eq. so the hadronic tensor 
+ iWr) reads 

Wr = -0(g°) I ^ Im[n(V" + n^^](g;p) (8) 

wr = -0(9°) / Cr 

from where we can see how is the basic object of our approach. Following the lines of Q, we should perform 
a many body expansion, where the relevant gauge-boson absorption modes would be systematically incorporated: 
absorption by one, two or even three nucleon mechanisms, real and virtual meson (tt, p, • • •) production, excitation 
of A or higher resonance degrees of freedom, etc. Some of these modes are depicted in Fig. [3 


III. QE CONTRIBUTION TO 


The virtual W'^ can be absorbed by one nucleon leading to the QE contribution of the nuclear response function. 
Such a contribution corresponds to a Iplh nuclear excitation (first of the diagrams depicted in Fig. |2J). We will work 
on a non-symmetric nuclear matter with different Fermi sea levels for protons than for neutrons. For the W^-pn 
vertex we consider the V — A current, and use PCAC and invariance under G-parity to relate the pseudoscalar form 
factor to the axial one and to discard a term of the form -|- in the axial sector, respectively. Invariance 

under time reversal guarantees that all form factors are real. Using isospin symmetry we can relate the vector form 
factors with the electromagnetic ones. 

With all of these ingredients is straightforward to evaluate the contribution to the W~^ self-energy of the first 
diagram of Fig. mwe finally get 


W>^'^{q\q) = - 


cos^ 9c 
2M2 


drr^2Q{q°) 


(Pp M M 

(27r)3 EpEp+ii 


(-tt) 


( 10 ) 


X 0(fc^ - Ip I)0(|p-bg I - kPp)S{q° + Ep- Ep+q)A’'f^{p,q)\pO^Ef; 


with the CC nucleon tensor obtained after taking some traces on the Dirac’s space. The cPp integrations above 
can be done analytically and all of them are determined by the imaginary part of the relativistic isospin asymmetric 
Lindhard function, U}i{q, kp, kp). Explicit expressions for Ur and are given in Q. 















4 



V 



V 





FIG. 3: Set of irreducible diagrams responsible for the polarization (RPA) effects in the Iplh contribution to the self-energy. 


Up to this point the treatment is fully relativistic. To account for RPA effects, we will use a nucleon-nucleon 
effective force, so for consistency we ought to use a non-relativistic Fermi gas. This is easily done by replacing the 
factors M/Ep and M/Ep^p in Eq. lUITI) by one. 

Pauli blocking, through the imaginary part of the Lindhard function, is the main nuclear effect included in the 
hadronic tensor of Eq. m- In the next sections we will study additional nuclear corrections to . 

A few words here on the low density theorem (LDT): when low nuclear density is supposed, the imaginary part of 
the Lindhard function can be approximated by a Dirac delta on energy (up to a constant factor) in such a way that 
the model reproduces the free space nucleon cross section. 


IV. NUCLEAR MODEL CORRECTIONS 
A. Nuclear Correlations 

When the electroweak interactions take place in nuclei the strengths of electroweak couplings may change from 
their free nucleon values due to the presence of strongly interacting nucleons; indeed, since the nuclear experiments 
on j3 decay in the early seventies |^, the quenching of axial current is a well established phenomenon. We follow here 
the many body approach of Q , and take into account the medium polarization effects in the Iplh contribution to the 
W~^ boson self-energy by substituting it with an RPA response as shown diagrammatically in Fig.|3For that purpose 
we use an effective ph-ph contact interaction 

v= Co {fo(p) + fo(p)nT2 + go(p)^1^2 + go(p)^1^2riT2} (11) 

of the Landau-Migdal type. The density dependent coefficients were determined [T^ from calculations of nuclear 
electric and magnetic moments, transition probabilities, and giant electric and magnetic multipole resonances. In the 
S = T = 1 channel (ctcttt operator) we use an interaction with explicit tt (longitudinal) and p (transverse) exchanges, 
which has been used for the renormalization of the pionic and pion related channels in different nuclear reactions at 
intermediate energies Hi Further effects such as short range correlations (SRC) are also taken into account. 

We also include A(1232) degrees of freedom in the nuclear medium which, given the spin-isospin quantum numbers 
of the A resonance, only modify the vector-isovector (S = T = 1) channel of the RPA response function. 

The V lines in Fig. 0 stand for the effective ph(Ah)-ph(Ah) interaction described so far. We should stress that 
this effective interaction is non-relativistic, and then for consistency we will neglect terms of order /NP) when 
summing up the RPA series. 
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FIG. 4; W'*' self-energy diagram obtained from the first diagram depicted in Fig. |21by dressing np the nucleon propagator of 
the particle state in the ph excitation. 


B. Energy Balance and Coulomb Distortion 

To ensure the correct energy balance in the reaction o for finite nuclei, the energy conserving Dirac delta function 
in Eq. m has to be modified by including the minimum excitation energy, Q = M{Az+i) — M{Az)i needed for 
the transition to the ground state of the final nucleus. The consideration of this energy gap is essential to obtain 
reasonable cross sections for low-energy neutrinos, see 0 - 

We also include a Coulomb self-energy Sc = 2k'°Vc{r) in the intermediate lepton propagator of the neutrino 
self-energy depicted in Fig.^where Vc{r) is the nucleus Coulomb potential produced by a charge distribution Pch{r). 
This way of taking into account the Coulomb effects has clear resemblances with what is called “modified effective 
momentum approximation” in 0 


C. ESI 

Once a ph excitation is produced by the virtual FT—boson, the outgoing nucleon can collide many times, thus 
inducing the emission of other nucleons. The result is a quenching of the QE peak respect to the simple ph excitation 
calculation and a spreading of the strength, or widening of the peak. A distorted wave approximation with an optical 
(complex) nucleon-nucleus potential would remove all these events. However, if we want to evaluate the inclusive 
l~) cross section these events should be kept and one must sum over all open final state channels. 

We will account for the Final State Interaction (ESI) by using nucleon propagators properly dressed with a realistic 
self-energy in the medium, which depends explicitly on the energy and the momentum [l3j. This self-energy has an 
imaginary part from the coupling to the 2p2h components, which is equivalent to the use of correlated wave functions, 
evaluated from realistic NN forces and incorporating the effects of the nucleon force in the nucleon pairs. Thus, we 
consider the many body diagram depicted in Fig. 

Once we have got a model for the nucleon self-energy E(p'',p;p), we can include in it a renormalized nucleon 
propagator Cpsiip'^p), that can be easily related to Sp and Sh, the particle and hole spectral functions and then 
through 


Im[/ FSiiq', kp) 


e(g°) 

47r^ 



dujSh{u!,p; p)Sp{q° 


+ uj,p + q\p) 


( 12 ) 


with the Lindhard function that we include in our formalism. 
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TABLE I: Flux averaged ^^C(i/e,e )X and )X cross sections. 


Theory 

KARMEN [IJ 

LSND [IJ 

LAMPF [16j 

CTe 0.14 

0.15 ±0.01 ±0.01 

0.15 ±0.01 ±0.01 

0.141 ±0.023 

Theory 

LSND’95 

LSND’97 

LSND’02 [17] 

11.9 

8.3 ± 0.7 ± 1.6 

11.2 ±0.3 ±1.8 

10.6 ±0.3 ±1.8 





FIG. 5: Predictions for the LSND experiment. See 0 for details. 


V. LOW ENERGY RESULTS 

We present in Figs. 0 and El and Table Q] our theoretical predictions and a comparison with the experimental 
measurements of the inclusive ii~)X and reactions near threshold. Pauli blocking and the use 

of the correct energy balance improve the results, but only once RPA and Coulomb effects are included a good 
description of data is achieved. 

Given the succes of the model at low energies we decided to further test it by calculating inclusive muon capture 
rates in nuclei throughout the Periodic Table. Results are given in Table im including the error in the theoretical 
predictions. Data were taken from [^, using a weighted average: F/ct^ = Yl,i l/cr,^. 


TABLE II: Experimental and theoretical total muon capture widths L, for different nuclei. See 0 for details. 



Pauli±(3 RPA 

Exp 



5.42 

3.21 

3.78 ± 0.03 

0.15 

16 0 

17.56 

10.41 

10.24 ±0.06 

-0.02 

ISq 

11.94 

7.77 

8.80 ±0.15 

0.12 

Na 

58.38 

35.03 

37.73 ±0.14 

0.07 

«Ca 

465.5 

257.9 

252.5 ±0.6 

-0.02 

^"Ca 

318 

189 

179 ±4 

-0.06 

As 

1148 

679 

609±4 

-0.11 

“2Cd 

1825 

1078 

1061±9 

-0.02 

208pb 

1939 

1310 

1311±8 

0.00 
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FIG. 6: Predictions for the )X reaction. See 0 for details. 


Despite the huge range of variation of the capture widths, the agreement to data is quite good for all studied nuclei, 
with discrepancies of about 15% at most. Furthermore, using LFG instead of a more refined model such as a shell 
model does not affect much the value of integrated observables such as total cross section or capture widths, see 
It is precisely for where we find the greatest discrepancy with experiment. Nevertheless, our model provides one 
of the best existing combined description of the inclusive muon capture in and the LSND measurement of the 
reaction /j,“)V near threshold. 


VI. INTERMEDIATE ENERGY RESULTS 

At intermediate energies the predictions of this model should become reliable, not only for integrated, but also for 
differential cross sections. We present results for incoming neutrino energies within the interval 150-400 (250-500) 
MeV for electron (muon) species. The use of relativistic kinematics for the nucleons leads to moderate reductions in 
the interval of 4-9% for both neutrino and antineutrino cross sections, at the energies considered. Such corrections 
do not depend significantly on the considered nucleus. 

In Fig.[7|the effects of RPA and Coulomb corrections are studied as a function of the incoming neutrino/antineutrino 
energy. The correction Rrpa+cb is defined as (ctrpa+cb — o’o)/croj where uq does not include RPA and Coulomb 
corrections, while ctrpa+cb includes these nuclear effects. FSI corrections are not taken into account in these cross 
sections. RPA correlations reduce the cross sections, and we see large effects, specially at lower energies. Nevertheless, 
for the highest energies considered (500 and 400 MeV for muon and electron neutrino reactions, respectively) we still 
find suppressions of about 20-30%. Coulomb distortion of the outgoing charged lepton enhances (reduces) the cross 
sections for neutrino (antineutrino) processes and its effects decrease with energy. For antineutrino reactions, the 
combined effect of RPA and Coulomb corrections have a moderated dependence on A and Z. At the high energy 
end the A—dependence becomes milder, since Coulomb distortion becomes less important. In the case of neutrinos, 
the increase of the cross section due to Coulomb cancels out partially with the RPA reduction. Finally, the existing 
differences between electron and muon neutrino/antineutrino plots are due to the different momenta of an electron 
and a muon with the same energy. 

In Figs. IHl and m FSI effects on differential cross sections are shown. As expected, FSI provides a broadening and 
a significant reduction of the strength of the QE peak. Finally, in Table cni we compile muon and electron neutrino 
and antineutrino inclusive QE integrated cross sections from oxygen. We present results for relativistic (REL) and 
non-relativistic nucleon kinematics and in this latter case, we present results with (FSI) and without FSI (NO-REL) 










FIG. 7: RPA and Coulomb (CB) corrections to electron and muon neutrino and antineutrino QE cross sections for different 
nuclei, as a function of the neutrino energy. 


effects. Though FSI changes importantly the shape of the differential cross sections, it plays a minor role when 
one considers total cross sections. When medium polarization effects are not considered, FSI provides significant 
reductions (13-29%) of the cross sections [^. However, when RPA corrections are included, the reductions becomes 
more moderate, always smaller than 7%; even there exist some cases where FSI enhances the cross sections. This 
can be easily understood by looking at Fig. where we show the differential cross section as a function of the energy 
transfer for = 375 MeV. There, we see that FSI increases the cross section for high energy transfer. But for nuclear 
excitation energies higher than those around the QE peak, the RPA corrections are certainly less important than in 
the peak region. Hence, the RPA suppression of the FSI distribution is significantly smaller than the RPA reduction 
of the distribution determined by the ordinary Lindhard function. 


VII. PREVIOUS RESULTS 

The same formalism presented here has been used in previous works studying real and virtual P| photon inclusive 
nuclear reactions. Excellent results both in the quasielastic and A excitation regions where obtained in these works. 
To describe the A peak and the ’’dip” regions, they included a high number of gauge boson absorption modes so they 
were able to study the reaction at higher nuclear excitation energies than those we have presented here. As can be 
seen for instance in Fig. HOI the agreement with experiment is excellent. Furthermore, inclusive processes of the type 
(e, e'A), {e,e'NN), (e,e'7r) were studied by means of a MonteCarlo simulation as presented in that make use of 
the nuclear and pion physics models of and 0 • It is because of the remarkable succes of this model that we expect 
our work to be highly reliable for CC reactions. 
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FIG. 8: Muon neutrino differential cross sections in calcium as a function of the lepton scattering angle (top) and of the 
momentum transfer (bottom). Top: Cross sections, without FSI and using relativistic kinematics for the nucleons. Crosses 
have been obtained without RPA and Coulomb effects, while the curves have been obtained with the full model (up to FSI 
effects). Bottom: Cross sections, obtained by using relativistic (REL) and non-relativistic nucleon kinematics results with 
(FSI) and without (NOREL) FSI effects. We also take into account RPA and Coulomb corrections (lower lines at the peak). 


VIII. CONCLUSION AND OUTLOOK 

We have presented here a many body approach to inclusive electroweak reactions in nuclei, at intermediate energies 
(nuclear excitation energies below 500 MeV). It systematically takes into account RPA, SRC, A(1232), FSI and MEC 
effects. The meson-nucleon and nucleon-nucleon dynamics of the approach have been successfully tested in former 
pionic reactions. 

It has been tested succesfully in; 

• Real and virtual photo-absorption and tt, V, NN, Nn electro and photoproduction processes in nuclei. 

• Charged current induced inclusive neutrino ^~)X cross sections at low energies and Inclusive Muon 

capture in Nuclei. 

Predictions for QE neutrino induced reactions in nuclei at intermediate energies of interest for future neutrino 
experiments have been presented. 

Our intention is to improve this approach by including contributions from resonance degrees of freedom and MEC 
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TABLE III: Muon and electron neutrino and antineutrino inclusive QE integrated cross sections from 




E, 


1 

O 






REL 

NOREL ESI 

REL 

NOREL ESI 

500 Pauli 

460.0 

497.0 

431.6 

155.8 

168.4 

149.9 

RPA 

375.5 

413.0 

389.8 

113.4 

126.8 

129.7 

375 Pauli 

334.6 

354.8 

292.2 

115.1 

122.6 

105.0 

RPA 

243.1 

263.9 

243.9 

79.8 

87.9 

87.5 

250 Pauli 

155.7 

162.2 

122.5 

63.4 

66.4 

52.8 

RPA 

94.9 

101.9 

93.6 

38.8 

42.1 

40.3 




1 

o 



0(/^e, 



REL 

NOREL ESI 

REL 

NOREL ESI 

310 Pauli 

281.4 

297.4 

240.6 

98.1 

104.0 

87.2 

RPA 

192.2 

209.0 

195.2 

65.9 

72.4 

73.0 

220 Pauli 

149.5 

156.2 

121.2 

60.7 

63.6 

51.0 

RPA 

90.1 

97.3 

92.8 

36.8 

40.0 

40.2 

130 Pauli 

37.0 

38.3 

28.8 

21.1 

21.9 

16.9 

RPA 

20.6 

22.3 

23.3 

10.9 

11.9 

12.8 



FIG. 9: Muon neutrino QE differential cross sections in as a function of the energy transfer. We show results for relativistic 
(long dashed line, ’REL’) and non-relativistic nucleon kinematics with (solid line, ’ESI’) and without (short dashed line, 
’NOREL’) ESI effects. We also show the effect of RPA and Coulomb corrections (lower lines at the peak). 


in the charged current reactions. We also want to extend this formalism to exclusive channels in neutral currents via 
a MonteCarlo simulation. 


IX. LIST OF SYMBOLS/NOMENCLATURE 

We have used h = c = 1 units for formulas all throughout this work, however results in tables and figures are 
presented in the following units unless otherwise noted. 

k = LAB lepton momenta, MeV Ep = Energy of p momentum lepton, MeV 

p = Nuclear matter density a = Cross section, 10“^°cm^ 

a = Spin Pauli matrices r = Isospin Pauli matrices 

fcp = Fermi momentum for neutrons kp = Fermi momentum for protons 
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FIG. 10: Double differential ^^C{e,e')X cross section 

q = Transferee! W momentum T = Muon capture widht, 10~^s“^ 

0(a:) = Step function 
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